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Summary 

The subsurface microbiology of an Athabasca oil 
sands reservoir in western Canada containing 
severely biodegraded oil was investigated by combin- 
ing 16S rRNA gene- and polar lipid-based analyses of 
reservoir formation water with geochemical analyses 
of the crude oil and formation water. Biomass was 
filtered from formation water, DNA was extracted using 
two different methods, and 16S rRNA gene fragments 
were amplified with several different primer pairs 
prior to cloning and sequencing or community finger- 
printing by denaturing gradient gel electrophoresis 
(DGGE). Similar results were obtained irrespective of 
the DNA extraction method or primers used. Archaeal 
libraries were dominated by Methanomicrobiales (410 
of 414 total sequences formed a dominant phylotype 
affiliated with a Methanoregula sp.), consistent with 
the proposed dominant role of C02-reducing methano- 
gens in crude oil biodegradation. In two bacterial 16S 
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rRNA clone libraries generated with different primer 
pairs, > 99% and 100% of the sequences were affiliated 
with Epsilonproteobacteria (n = 382 and 72 total 
clones respectively). This massive dominance of 
Epsilonproteobacteria sequences was again obtained 
in a third library (99% of sequences; n = 96 clones) 
using a third universal bacterial primer pair (inosine- 
341 f and 1492r). Sequencing of bands from DGGE 
profiles and intact polar lipid analyses were in 
accordance with the bacterial clone library results. 
Epsilonproteobacterial OTUs were affiliated with Sul- 
furicurvum, Arcobacter and Sulfurosplrillum spp. 
detected in other oil field habitats. The dominant 
organism revealed by the bacterial libraries (87% of all 
sequences) is a close relative of Sulfuricurvum 
kujiense - an organism capable of oxidizing reduced 
sulfur compounds in crude oil. Geochemical analysis 
of organic extracts from bitumen at different reservoir 
depths down to the oil water transition zone of these 
oil sands indicated active biodegradation of diben- 
zothiophenes, and stable sulfur isotope ratios for 
elemental sulfur and sulfate in formation waters were 
indicative of anaerobic oxidation of sulfur com- 
pounds. Microbial desulfurization of crude oil may be 
an important metabolism for Epsilonproteobacteria 
indigenous to oil reservoirs with elevated sulfur 
content and may explain their prevalence in formation 
waters from highly biodegraded petroleum systems. 

Introduction 

The global inventory of petroleum reserves is dominated 
by heavy oil reservoirs that represent a legacy of anaero- 
bic microbial communities that have degraded hydrocar- 
bons in situ over geological timescales (Head et al., 2003; 
2010; Larter etal., 2006). These biologically altered fossil 
fuels represent an enormous energy resource, yet their 
production is complicated by detrimental factors such as 
increased oil viscosity, acidity and sulfur content. Extreme 
cases of biodegradation in 'unconventional' super-heavy 
oil reservoirs have led to the development of novel extrac- 
tion techniques involving mining or steam injection for 
producing this petroleum. A prominent example is the 
Athabasca oil sands (Fig. 1), where increased production 
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excavation area 



Fig. 1. Geographic location of the Athabasca oil sands. Three formation water samples were taken from around the Muskeg River open pit 
mine site where unconventional heavy oil recovery operations take place. The approximate location of the Pelican Lake oil reservoir where 
conventional recovery is practised (Grabowski etal., 2005a,b) is also indicated (P). 



in recent years underscores Canada having the second 
largest proven domestic crude oil reserves worldwide (c. 
170 billion barrels; ERCB, 2009). The many challenges 
associated with producing heavy oil or bitumen (highly 
viscous oil) translate into added costs and increasing 
environmental impacts, and highlight the ongoing need for 
innovation in the oil sands sector. Biotechnologies based 
on understanding and manipulating subsurface reservoir 
microbial communities have been proposed as a promis- 
ing route towards more sustainable production of energy 
from fossil fuels including oil sands (Grigoryan and Voor- 
douw, 2008; Jones etal., 2008; Youssef etal., 2009). 
These strategies are underpinned by a view of petroleum 
reservoirs as important habitats within the 'deep bio- 
sphere'. Investigating reservoir habitats is often enabled 
by sampling produced water, which has been conducted 
for several oil fields including some in western Canada 
(Voordouw etal., 1996; Grabowski etal., 2005a). Similar 
studies from unconventional oil sands reservoirs have not 
been reported to date. 

The occurrence of fermentative bacteria and methano- 
genic archaea in several oil fields (Jeanthon etal., 2005; 
Ollivier and Cayol, 2005) befits current models of anaero- 



bic hydrocarbon degradation in subsurface petroleum res- 
ervoirs. Hydrocarbon geochemistry and gas isotope data 
have highlighted crude oil hydrocarbon degradation domi- 
nated by CO2 reduction to methane as an important 
mechanism (Jones etal., 2008; Gray etal., 2009; 2010). 
Experimental microcosms inoculated with surface sedi- 
ments show that anaerobic degradation of hydrocarbons 
and crude oil can be catalysed by consortia dominated 
by hydrogenotrophic (C02-reducing) methanogens (e.g. 
order Methanomicrobiales) and fermentative syntrophic 
Deltaproteobacteria (e.g. family Syntrophaceae) (Zengler 
etal., 1999; Jones etal., 2008). This hypothesis is con- 
sistent with the majority of the methanogenic archaea 
recovered so far from oil fields being hydrogenotrophic 
CO2 reducers (Jeanthon etal., 2005; Nazina etal., 2006; 
Gray etal., 2009). Interestingly, bacterial community 
analyses of reservoir samples have not revealed a clear 
prominence of Syntrophaceae to mirror the widespread 
occurrence of Methanomicrobiales. The absence, or low 
abundance of these putative syntrophs in oil reservoir 
communities may indicate that other bacterial groups 
might fill this niche in subsurface petroleum ecosystems. 
An alternative explanation could be the relatively few 
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investigations of indigenous microbial communities from 
actively biodegrading oil reservoirs where syntrophs are 
expected to be important in situ. Grabowski and col- 
leagues (2005a, b) examined formation waters (i.e. 
co-produced with crude oil during primary recovery opera- 
tions prior to secondary water flooding) from the Pelican 
Lake oil field - a 400 m deep conventional heavy oil 
reservoir in west Athabasca (Fig. 1). Fatty acid-degrading 
Syntrophus spp. could be enriched from this formation 
water after several transfers (Grabowski etal., 2005b), 
but were not detected in a bacterial 1 6S rRNA gene clone 
library constructed from DNA extracted from formation 
water (Grabowski etal., 2005a). This clone library con- 
sisted of 1 51 sequences from a single taxon affiliated with 
the genus Arcobacter within the Epsilonproteobacteria 
(Grabowski etal., 2005a). 

Even though Epsilonproteobacteria have been 
detected in several oil fields using a variety of methods 
(e.g. Voordouw etal., 1996; Telang etal., 1997; Gevertz 
etal., 2000; Grabowski etal., 2005a; Hubert and Voor- 
douw, 2007; Sette etal., 2007; Gittel etal., 2009; Pham 
et al., 2009), the unexpected predominance of sequences 
from Epsilonproteobacteria in the Pelican Lake clone 
library has led to the suggestion that this may represent 
an artefact arising from selective amplification of epsilon- 
proteobacterial 16S rRNA gene sequences by certain 
PGR primers (Grabowski etal., 2005a; Pham etal., 2009; 
Head etal., 2010). Such primer selectivity was shown to 
explain similar results obtained by Watanabe and col- 
leagues (2000) who reported strong dominance of Epsi- 
lonproteobacteria in a similar environment (91% of clones 
derived from groundwater in an underground crude oil 
storage cavity), but later reported that using different 
PGR primers resulted in recovery of a greater diversity 
of sequences and a smaller proportion affiliated with 
Epsilonproteobacteria (no more than 16% of clones; 
Watanabe etal., 2002). 

Here we present a microbiological characterization of 
formation waters from an unconventional shallow oil 
sands reservoir (0-80 m) from the northern Athabasca 
region in western Canada where the oil is severely 
degraded and immobile. Geochemical, microbiological, 
and lipid- and DNA-based molecular analyses of crude oil 
and formation water were performed, including PCR- 
based screening of 16S rRNA gene sequences from resi- 
dent Archaea and Bacteria by employing different DNA 
extraction and PGR protocols to circumvent selectivity 
imposed by using only a single approach. The results 
reveal a strong dominance of Epsilonproteobacteria in the 
formation water samples from this oil sands reservoir. 
Coupling this observation with geochemical data showing 
that dibenzothiophenes are extensively biodegraded 
throughout the oil column and particularly at the oil water 
transition zone, we suggest that the predominance of 



Epsilonproteobacteria may be explained by the ability of 
these bacteria to use organic sulfur compounds in crude 
oil as an electron donor and energy source (Kodama and 
Watanabe, 2003). 

Results 

Site description and formation water geochemistry 

Samples were obtained from the IVIuskeg River mine, 
located approximately 75 km north of the city of Fort 
IVlcMurray, Alberta (Fig. 1). In this oil sands reservoir, for- 
mation water occurs as a basal aquifer up to 20 m in 
thickness underlying a 50-80 m thick layer of oil- 
saturated sands. Formation waters are therefore not in 
contact with surface waters. Muskeg River oil sands are 
produced by clearing several metres of topsoil (overbur- 
den) to allow large-scale truck and shovel operation to 
excavate and transport oil sands to nearby facilities that 
separate the bitumen from the sand. This unconventional 
surface mining approach requires advanced dewatering 
of the oil sands to reduce formation pressure and prevent 
fractures and flooding as the open pit excavation pro- 
ceeds deeper into the reservoir. Therefore formation 
waters are discharged at dewatering wells located 500- 
2500 m ahead of the advancing production area. Dis- 
charged formation waters should thus derive from 
biologically active oil water transition zones and associ- 
ated aquifers in the subsurface (Head etal., 2003; Larter 
etal., 2003) and represent pristine reservoir samples 
unaffected by any prior water or chemical injections that 
sometimes characterize produced water samples from 
deeper conventional oil fields (IVIagot, 2005). Waters dis- 
charged from six wells around the mine excavation area 
had a combined flow rate of 250 m^ h"\ and three of these 
wells were sampled for geochemical and microbiological 
analyses (Fig. 1). The formation water geochemistry is 
summarized in Table 1, and indicates that the three 
samples are broadly similar. Similar microbial community 
compositions in the three water samples were confirmed 
by denaturing gradient gel electrophoresis (DGGE) (not 
shown) and intact polar lipid (IPL) analyses (Oldenburg 
et al., 2009). Further community analyses focused on for- 
mation water from wellhead #1 . 

Formation water microbial community composition 

Two archaeal clone libraries (one prepared in Newcastle 
and the other in Calgary and Rockville) and three bacterial 
clone libraries (two prepared in Newcastle and one in 
Calgary and Rockville) were constructed from DNA that 
was extracted from the biomass that had been concen- 
trated by filtering 1 2 I of formation water from wellhead #1 . 
The diversity in all libraries was low with just seven 
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Table 1. Geochemistry of oil sands formation water samples. 



Parameter 




Wellhead #1 


Wellhead #4 


Wellhead #6 




Basal aquifer thicl<ness 




18.6 m 


17.6 m 


23.0 m 




Cumulative discharge at the time of sampling 227 000 m^ 


373 000 m^ 


387 000 m= 




Discharge flow rate (flow continuity) 


10-20 m" h-' 


20-30 m^ h-' 


60-70 m" h-' 








(continuous) 


(semi-continuous) 


(semi-continuous) 


Electrical conductivity mS cm"' 


2.36 


5.22 


4.20 




pH 




7.20 


7.38 


7.43 




Alkalinity (mg M HCOa") 




959.3 


1512.0 


1387.0 




Intact polar lipids detected 




Phosphatidylethanolamine; Phosphatidylethanolamine; Phosphatidylethanolamine; 






phosphatidylglycerol 


phosphatidylglycerol 


phosphatidylglycerol 


S"S elemental S (%o) 




ND 


26.7 


26.7 




S"S SO/- (%=) 




22.8 


23.6 


23.4 




0U4 (.mg 1 ) 




15.37 


32.50 


23.03 




NO3- 




ND 


ND 


ND 




po/- 




ND 


ND 


ND 




Mn^* 




0.27 


0.04 


0.02 








0.50 


0.15 


0.10 




Na* 




454.00 


1276.00 


1033.10 




ci- 




351.22 


1226.11 


877.24 




Br 




5.71 


9.68 


8.60 




F- 




3.23 


3.60 


3.72 




Si* 




9.24 


4.33 


5.48 




Ba^+ 




0.81 


0.16 


0.22 








2.04 


1.84 


1.63 




ur 




0.22 


0.25 


0.22 








77.07 


51.11 


52.78 








43.42 


31.84 


30.70 




K* 




18.10 


22.10 


19.40 




archaeal and five bacterial phylotypes being recovered 


munity compositions were obtained for Archaea and 


from tlie five libraries 


(Table 2; operational taxonomic 


Bacteria in the oil sands formation water. The Newcastle 


units (OTUs) were defined at 97% 


sequence identity). 


archaeal library (91 clones) consisted of just a single 


Clone library results 


indicate that 


regardless of DNA 


phylotype affiliated with 


Methanoregula (Table 2). This 


extraction procedure or PGR primers used, similar com- 


phylotype was also dominant (89% of clones) in the larger 


Table 2. Archaeal and bacterial 16S rRNA gene clone library results. 








Library; target group 


No. of clones 










(16S primers) 


sequenced 


Lineages of 


% abundance in 


Representative 


Accession 


(DNA extraction) 


in total 


phylotypes detected 


library (No. of clones) 


type sequence 


No. 


A1 Archaea 


323 


Methanoregula 


88.9% (287) 


TS1A275 


JF789587 


(arch8f/arch1492r) 




Methanospirillum 


9.3% (30) 


TS1A121 


JF789588 


(Calgary) 




Methanomicrobiales 


0.6% (2) 


TS1A142 


JF789589 






Methanosarcina 


0.3% (1) 


TS1A083 


JF789590 






Desulfurococcales 


0.3% (1) 


TS1A038 


JF789591 






Thermoplasmatales 1 


0.3% (1) 


TS1A251 


JF789592 






Thermoplasmatales 2 


0.3% (1) 


TS1A042 


JF789593 


A2 Archaea 


91 


Mettianoregula 


100% (91) 


(TS1A275)= 


(JF789587) 


(arch46f/arch1017r) 












(Newcastle) 












B1 Bacteria (9f/1545r) 


382 


Sulfuricurvum 


95.6% (369) 


TS1B301 


JF789594 


(Calgary) 




Arcobacter 


2.9% (11) 


TS1B220 


JF789595 






Suifurospirillum 1 


0.3% (1) 


TS1B252 


JF789596 






Peiobacter 


0.3% (1) 


TS1B322 


JF789597 


B2 Bacteria (?,\l-\5A2r) 


72 


Suifuricurvum 


63.8% (46) 


(TS1B301) 


(JF789594) 


(Newcastle) 




Arcobacter 


31 .9% (23) 


(TS1 B220) 


(JF789595) 






Suifurospirillum 2 


4.2% (3) 


NCL08_D6E05 


JF789598 


B3 Bacteria 


96 


Sulfuricurvum 


63.5% (61) 


(TS1B301) 


(JF789594) 


(ino-341f/1492r) 




Arcobacter 


29.2% (28) 


(TS1 B220) 


(JF789595) 


(Newcastle) 




Suifurospirillum 2 


6.3% (6) 


(NCL08_D6E05) 


(JF789598) 






Peiobacter 


1.0% (1) 


(TS1 B322) 


(JF789597) 



a. Designations in parentheses indicate the type sequence is from one of the other libraries. 
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J. Craig Venter Institute (JCVI) library (n = 323 clones 
analysed), which contained additional phylotypes related 
to l\/letiianospirilium (9% of clones), unclassified Metlia- 
nomicrobiales (0.6%), Metlianosarcina (0.3%), Desulfuro- 
coccales (0.3%) and Tliermopiasmatales (0.3% for each 
of two distinct OTUs) (Table 2). 

Common to all three bacterial libraries were two 
phylotypes related to Sulfuricurvum and Arcobacter, 
with Sulfuricurvum always being the most abundant 
group (Table 2). Arcobacter comprised a much higher 
proportion of clones in Newcastle libraries, which were 
constructed from DNA that was extracted using a bead 
beating method (30% Arcobacter clones in Newcastle 
libraries as opposed to 3% in the JCVI library; Table 2). 
Both of the Newcastle libraries revealed a distinct 
Sulfurospirillum phylotype (4% and 6% of clones) 
that was not detected in the larger JCVI library, and the 
JCVI library contained a distinct Sulfurospirillum phylo- 
type (a singleton) that was not detected in the Newcastle 
libraries. The only non-epsilonproteobacterial clones 
were singletons from the two larger libraries that com- 
prised a fourth phylotype related to Pelobacter in the 
Deltaproteobacteria. 

Only two bacterial clone libraries were constructed ini- 
tially, using primers 9f/1545r and 8f/1542r (Table S1), 
resulting in 99.7% and 1 00% of the cloned 1 6S rRNA gene 
sequences being affiliated with Epsilonproteobacteria 
(Table 2). This was unexpected, and similar to the 91% 
dominance of Epsilonproteobacteria in an underground 
storage cavity reported by Watanabe and colleagues 
(2000) using a similar primer pair (8f/1546r). Watanabe 
and colleagues (2002) subsequently constructed an addi- 
tional library using an inosine-substituted 341 f primer 
together with 1492r, which resulted in only 16% of clones 
affiliated with Epsilonproteobacteria (also related to Sulfu- 
ricurvum, Arcobacter and Sulfurospirillum). However, 
when the primers designed by Watanabe and colleagues 
(2002) were used to construct a third bacterial clone 
library from the oil sands formation water, again 99% of 
sequences were affiliated with Epsilonproteobacteria 
(Table 2). 

The bacterial community composition determined 
by the three clone libraries was confirmed using a 
nested DGGE strategy whereby different universal bac- 
terial primer pairs 8f/1542r or inosine-341f/1492r were 
used to amplify 16S rRNA genes from the formation 
water DNA sample, with resulting PCR products used 
as a template for DGGE-PCR using 341f-GC/534r 
(Table SI). Figure 2 shows that both methods resulted in 
the same banding pattern. All bands that were excised 
and sequenced corresponded to Epsilonproteobacteria, 
having 97-100% sequence identity to the dominant 
Sulfuricurvum phylotype from the clone libraries (data 
not shown). 



i-341f 8f 
M 1492r M 1542r M 




Fig. 2. Denaturing gradient gel electrophoresis of bacterial 16S 
rRNA gene fragments following two-step nested PCR. Amplification 
with different universal bacterial primers (indicated above gel lanes; 
see also Table St) was followed by a second round of amplification 
using DGGE primers to create ampllcons 233 base pairs in length 
(including a 40-base-pair GC clamp) for DGGE, as described In the 
text. Arrows indicate bands from lanes 2-3 and 5-6 that were 
excised and sequenced. All sequences were closely related to the 
Sulfuricurvum OTU detected in the bacterial 16S rRNA gene clone 
libraries. 



Archaea detected in oil sands formation waters 

Putative hydrogenotrophic methanogens from the order 
Methanomicrobiales accounted for 410 out of 414 
archaeal sequences retrieved from the oil sands forma- 
tion water, comprising two dominant phylotypes related 
to Methanoregula and Methanospirillum spp. (Table 2; 
Fig. 3). Closely related Methanoregula spp. (> 98%; 
Fig. 3) have been detected in lake bottom sinkholes 
(e.g. Nold etal., 2010), a TCE-contaminated aquifer 
(Macbeth etal., 2004), a toluene-degrading enrich- 
ment (Ficker etal., 1999) and peatlands (e.g. Cadillo- 
Quiroz etal., 2008). One such peat bog was the source 
for the isolation of the type species Methanoregula 
boonei (Brauer etal., 2006), an H2- and C02-utilizing 
acidophilic methanogen that does not use other sub- 
strates (e.g. acetate, formate or trimethylamine; Brauer 
etal., 2011) and is the closest cultured relative to the 
dominant oil sands phylotype (Fig. 3). The second 
most abundant phylotype has the hydrogenotrophic 
Methanospirillum hungatei as its closest cultured rela- 
tive. Closely related Methanospirillum spp. have also 
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84j- Lake Huron sinkhole clone F18 (Nold et al., 2010) EU910620 
99P — TCE-contaminated groundwater clone TANA6 (Macbeth et al., 2004) AY667274 

L oil sands Methanoregula (92% of clones) JF789587 

I — municipal wastewater clone 57-1 (unpublished) AF424773 
'^_|- low H2 co-culture clone SP-ProM-A (Sakai et al., 2009) AB236066 

Methanoregula formicicum strain SMSP (Yashiro et al., 2011) AB479390 

peatland clone MHLsu47_10H (Cadillo-Quiroz et al., 2008) EU155920 

contaminated aquifer clone WCHD3-07 (Dojka et al., 1998) AF050617 

— toluene-degrading methanogenic enrichment clone M1 (Picker et al., 1999) AF423187 
Methanoregula boonel strain 6A8 (Brauer et al., 2006) DQ282124 

— hexadecane-degrading methanogenic enrichment clone A3 (Zengler et al., 1 999) AJ1 33793 
Pelican Lake oil field formation water clone PL-10E3 (Grabowski et al., 2005a) AY570644 

— Methanoculleus receptacull (T) isolated from an offshore oil field (Cheng et al., 2008) DQ787476 

— Pelican Lake oil field methanogenic enrichment clone FA-PA1 (Grabowski et al. 2005b) AY586393 

— Methanocalculus halotolerans isolated from an onshore oil field (Ollivier et al., 1998) AF033672 
Pelican Lake oil field formation water clone PL-10F8 (Grabowski et al., 2005a) AY570674 

iQQl Methanocorpusculum labreanum isolated from a tar pit (Zhao et al., 1989) AY260436 




Pelican Lake oil field formation water clone PL-9G3 (Grabowski et al., 2005a) AY570674 
Methanosplhllum hungatei (T) M60880 
hexadecane-degrading methanogenic enrichment clone A2 (Zengler et al., 1999) AJ133792 
- lake sediment clone Pav-sed-202 (unpublished) GU1351463 

■ peatland clone MH1492_11H (Cadillo-Quiroz et al., 2008) EU155983 

I oil sands Methanospirillum (7% of clones) JF789588 

Lake Huron sinkhole clone F5 (Nold et al., 2010) EU910621 
r TCE-contaminated groundwater clone TANAS (Macbeth et al., 2004) AY667273 
'°°L Pelican Lake sulfate-reducing enrichment clone PL-21A5 (Grabowski et al., 2005a) AY570675 

p oil sands Methanomicrobiales (0.5% of clones) JF789589 

Laptev Sea permafrost sediment clone PASLSS0.5m_25 (unpublished) FJ982770 
.mo Pelican Lake oil field formation water clone PL-10H6 (Grabowski et al., 2005a) AY570652 

' Pelican Lake oil field methanogenic enrichment clone FA-PA5 (Grabowski et al., 2005b) AY586394 

hexadecane-degrading methanogenic enrichment clone Al (Zengler et al., 1999) AJ 133791 

Pelican Lake oil field enrichment clone PL-37A1 (Grabowski et al., 2005a) AY570358 

oil sands Methanosarcina (0.2% of clones) JF789590 

^ Methanosarcina lacustrls (Simankova et al., 2003) AY260431 

oil sands Desulfurococcales (0.2% of clones) JF789591 

limestone sinkhole clone LPBB80 (unpublished) FJ901719 

100^ oil sands Thermoplasmatales 1 (0.2% of clones) JF789592 

' Antarctic cold seep clone ANTXXIII_706-4_Arch282 (unpublished) FN429784 

oil sands Thermoplasmatales 2 (0.2% of clones) JF789593 

Lake Huron sinkhole clone HI 3 (Nold et al., 2010) EU910628 

Fig. 3. Neighbour joining phylogenetic tree indicating taxonomic affiliations of cloned archaeal 16S rRNA gene sequences from Athabasca oil 
sands formation waters. For each of the archaeal OTUs the most closely related sequences from GenBank and other related sequences of 
interest (e.g. from other petroleum reservoirs or methanogenic hydrocarbon-degrading systems) are indicated. Bootstrap values next to 
branching nodes are based on 100 resamplings. 



been detected in sinkholes, aquifers and peatlands 
(see Fig. 3), and also include a Methanospirillum sp. 
enriched fronn the Pelican Lake oil reservoir c. 200 km to 
the west of these oil sands (Grabowski etal., 2005a; 
Figs 1 and 3). One out of the 414 sequences retrieved 
was affiliated with the Metiianosarclnales (Table 2; 
Fig. 3) and is closely related to a Metiianosarcina sp. 
detected In an archaeal 16S rRNA gene clone library 
from Pelican Lake formation water (without prior selec- 
tive enrichment). Three more singletons from the 



archaeal libraries were affiliated with non-methanogen 
sequences retrieved from cold seeps and sinkholes 
(Table 2; Fig. 3). 

Bacteria detected in oil sands formation waters 

Bacterial sequences affiliated with Sulfuricurvum were 
most closely related to environmental clones derived 
from meromictic lake sediments (Nelson etal., 2007), 
chemoclines (unpublished GenBank entry Accession No. 
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GQ390209; see Fig. 4), and sulfidic caves and springs 
(Porter and Engel, 2008) (Fig. 4). The only cultivated Sui- 
furicurvum reported so far are strains of Sulfuricurvum 
l<ujiense isolated from an underground crude oil storage 
cavity at Kuji, Japan (Kodama and Watanabe, 2003; 
2004), which have high 16S rRNA sequence identity to 
the oil sands formation water Sulfuricurvum sp. 16S rRNA 
sequence (98%; Fig. 4). Sulfuricurvum l<ujiense strain 
YK-1 is considered an obligate chemolithotroph (Kodama 
and Watanabe, 2004) that can grow on crude oil by 
coupling nitrate reduction to the oxidation of reduced, 
presumably organic, sulfur compounds in petroleum 
(Kodama and Watanabe, 2003). Sulfuricurvum kujiense 
cannot use hydrocarbons directly as a carbon and energy 
source but is able to oxidize sulfide, elemental sulfur, 
thiosulfate and hydrogen (Kodama and Watanabe, 2003; 
2004). 

Arcobacter s'p'p. are found in various habitats (Gevertz 
etal., 2000; Campbell e?a/., 2006; Miller etal, 2007; 
Fedorovich etal., 2009; Webster etal., 2010), and those 
closely related to the oil sands phylotype come from a 
range of environments (Fig. 4), including the Kuji oil cavity 
(Watanabe et al., 2000). The most closely related cultured 
Arcobacter (strain R-28314; Fig. 4) was isolated from a 
denitrifying wastewater treatment plant (Heylen etal., 
2006). An Arcobacter XhaX made up 100% of the bacterial 
16S rRNA gene clone library in Pelican Lake formation 
waters (Grabowski etal., 2005a) shares 92% 16S rRNA 
sequence identity with the Arcobacter detected in the oil 
sands formation waters studied here (Fig. 4). 

The two Sulfurospirillum phylotypes discovered in this 
study share 94% 16S rRNA sequence identity with each 
other (Fig. 4). Sulfurospirillum phylotype 2 has among its 
closest relatives an organism enriched from the Pelican 
Lake reservoir formation water in medium selective for 
fermentative bacteria (Grabowski etal., 2005a). Other 
relatives include uncultured bacteria from a TCE- 
contaminated subsurface aquifer (Macbeth etal., 2004) 
and soil associated with a Himalayan glacier (unpublished 
GenBank entries, e.g. Accession No. EU809872; Fig. 4). 
The closest cultivated relatives of both Sulfurospirillum 
phylotypes are Sulfurospirillum multivorans and S. halo- 
respirans, which are metabolically versatile facultative 
chemolithotrophs (Luijten etal., 2003; Fig. 4). Sulfurospir- 
illum phylotype 1 is most closely related (94%) to an 
uncultured lacustrine organism (unpublished GenBank 
entry Accession No. FJ612310; Fig. 4). Sulfurospirillum 
spp. have been detected in different oil field habitats, 
including the Kuji oil cavity. Pelican Lake and other fields 
in western Canada (Grabowski etal., 2005a; Hubert and 
Voordouw, 2007; Kodama etal., 2007; Fig. 4). 

The Pe/ofeacfer sequence detected in the larger bacte- 
rial clone libraries (Table 2) is most closely related to 
uncultured Deltaproteobacteria from different environ- 



ments including the sulfidic zone of an anoxic aquifer 
contaminated with heavy oil (Winderl etal., 2008). The 
most closely related described strain is Pelobacter propl- 
onicus strain DSM 2379 (Fig. 4), which ferments alcohols 
to organic acids (Schink etal., 1987). 

Intact polar lipids 

Microbial community analysis of the formation water 
samples also included isolation and characterization of 
IPLs from biomass that was filtered from the three forma- 
tion water samples. Phosphatidylethanolamines (PE) and 
phosphatidylglycerols (PG) were detected in all three 
samples (Table 1 ; for further details see Oldenburg et al., 
2009), indicating that formation waters retrieved from 
distant dewatering wells (kilometres apart; Fig. 1) may 
harbour similar bacterial communities. IPL structures can 
be used to make taxonomic inferences (Sturt et al., 2004), 
and while this approach lacks the phylogenetic resolution 
of 16S rRNA-based comparisons, bulk IPL analyses are 
not prone to issues such as primer specificity that can 
result in selective amplification of 16S rRNA genes from 
certain taxa. PE and PG are not highly diagnostic, as they 
are characteristic of many bacteria including Proteobac- 
teria (Dowhan, 1997). In a recent study PE and PG were 
reported to be the most abundant IPLs in the anoxic zone 
of the Black Sea water column (Schubotz etal., 2009) 
where nucleic acid-based cloning and fluorescence in situ 
hybridization analyses both revealed dominance of Epsi- 
lonproteobacteria (Vetriani etal, 2003; Lin etal., 2006), 
and high rates of sulfide oxidation are coupled to nitrate 
reduction (Jorgensen etal., 1991; Wakeham etal., 2007). 
We did not detect IPLs characteristic of archaea in oil 
sands formation waters, suggesting that bacteria are 
more abundant than archaea in these samples. 

Sulfur biogeochemistry in oil sands reservoir 
formation waters 

During sampling of the reservoir formation waters, a dis- 
tinct smell of hydrogen sulfide was noticed. The hydrogen 
sulfide concentration in the formation water samples was 
not measured during field sampling, but later chemical 
analysis detected up to 0.3 mM sulfate in the water 
samples (Table 1). Sulfur isotope analysis revealed mean 
6^''S values for elemental sulfur and sulfate of +26.7%o and 
+23.3%(j, respectively (Table 1 ), suggesting sulfate gen- 
eration via microbial oxidation of elemental sulfur and 
possibly other reduced sulfur compounds. This 4%o differ- 
ence in d^^S between sulfur and sulfate is consistent with 
isotope fractionation effects observed in a pure culture of 
sulfide-oxidizing Epsilonproteobacteria {Sulfurimonas sp. 
strain CVO) from a western Canadian oil field, grown 
under nitrate-reducing conditions (Hubert etal., 2009). 
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various oil field Sulfurospihtlum spp. (e.g., Grabowski et al., 2005a; Hubert & Voordouw, 2007) 

Dongping Lake done DP10.1.1 (unpublished) FJ612310 

oil sands Sulfurospirillum 1 (0.2% of clones) JF789596 

SulfurospiriHum cavolei isolated from Kuji oil storage cavity (Kodama et al., 2007) AB246781 
78^" Pelican Lake fermentative enrichment clones (Grabowski et al., 2005a) e.g., clone PL-26B3; AY570661 
SulfurospiriHum multivorans (T) (Scholz-Muramatsu et al., 1995) X82931 
. I — Sulfurospirillum tialorespirans (Luijten et al., 2003) NR_028771 
- Kuji oil storage cavity clone 1014 (Watanabe et al., 2000) AB030587 
glacier soil clones (unpublished) e.g., clone KS-318; EU809872 
acier soil clone KF-7 (unpublished) EF445480 
TCE-contaminated groundwater clone TAN6B (Macbeth et al., 2004) AY667251 
oil sands Sulfurospirillum 2 (1.6% of clones) JF789598 

methane seep sediment enrichment clone Mn3b-B13 (Beal et al., 2009) FJ264592 

Benzene-degrading enrichment clones (Herrmann et al., 2010) e.g., clone ZZ-S4D12; EF613424 
Arcobacter siram FWKOB isolated from a Canadian oil field (Gevertz et al., 2000) AF144693 
Kuji oil storage cavity clones (Watanabe et al., 2000) e.g., clone 1006; AB030597 
estuarine sediment clone 13CSRZ-B1 (Webster et al., 2010) FN424330 
western Canadian oil field clone ARC04 (Voordouw et al., 1996) U46505 

Pelican Lake oil field formation water clones (Grabowski et al., 2005a) e.g., clone PL-8B1; AY570593 
Alaska North Slope oil field clone D004025C08 (Pham et al., 2009) EU721819 
Arcobacter strain ED-1 isolated from a microbial fuel cell (Fedorovich et al., 2009) FJ968634 
Brazilian offshore oil field clone bo13C09 (Sette et al., 2007) AY862492 
North Sea oil f eld clone □anBac25 (Gittel et al., 2009) FN356234 
Black Sea sediment clone A3b2 (Thamdrup et al., 2000) AJ271655 

Arcobacter strain R-28314 isolated from a denitrifying wastewater (Heylen et al., 2006) AM084114 
i2p wastewater treatment plant biofilm clone 156 (Fernandez et al., 2008) AY692044 
~L Alaska North Slope oil field clone D004022A07 (Pham et al., 2009) EU721785 
oil sands /trcobacfer (11% of clones) JF789595 
vegetable harvest processing clone ATB-LH-6148 (unpublished) FJ535178 
- human liver carcinoma clone (unpublished) AF246706 
Kuji oil storage cavity clone 1053 (Watanabe et al., 2000) AB030593 
Sulfurimonas strain CVO isolated from a Canadian oil field (Gevertz et al., 2000) U46506 
estuarine sediment clone 13CSRZ-B8 (Webster et al., 2010) FN424332 
gold mine water microbial mat clone ZSB-D4-6 (unpublished) GU205603 
Sulfuricun/um kujiense strain YK-1 isolated from Kuji oil storage cavity (Kodama & Watanabe, 2003) AB053951 
- Frassassi cave spring clone FC5A_EPS2_cl3 (Porter & Engel, 2008) DQ295663 

J Frassassi cave water clone PC06110_B17 (Macalady et al., 2008) EU101169 

T-d Kuji oil storage cavity clones (Watanabe et al., 2000) e.g., clone 1005; AB030600 
< oil sands Sulfuricurvum (87% of clones) JF789594 
meromictic lake water clone SLB117 (Nelson et al., 2007) DQ787677 
- meromictic lake chemocline clone eub62B6 (unpublished) GQ390209 
aphotic cave sulfidic springs clone LKC3_127B.26 (Engel et al., 2004) AY510177 

Clones from underground caves and springs (Porter & Engel, 2008) e.g., clone CC1_cL51 ; DQ295579 
Smithella propionica strain LYP (Liu et al., 1999) AF1 26282 
hexadecane-degrading methanogenic enrichment clone B1 (Zengler et al., 1999) AJ133794 

Pelican Lake oil field methanogenic enrichment clone (Grabowski et al., 2005b) AY586399 

^ North Sea oil field formation water clone TCB131x (Dahle et al., 2008) DQ647155 

' Dutch onshore oil field clone c95Bk26b (van der Kraan et al., 2010) FJ941547 

Alaska North Slope oil field clone D004025D03 (Pham et al., 2009) EU721785 
Pelobacter propionicus DSM 2379 (Evers et al., 1993) X70954 
heavy oil contaminated aquifer clone D12_07 (Winderl et al., 2008) EU266817 
freshwater seep clone IS-83 (Bruun et al., 2010) GQ339169 
oil sands Pelobacter (0.2% of clones) JF789597 

C"^ pharmaceutical wastewater clone C11 (unpublished) EU234232 



meromictic lake chemocline clone HH285 (Halm et al., 2009) FJ502274 

Fig. 4. Neighbour joining phylogenetic tree indicating taxonomic affiliations of cloned bacterial 16S rRNA gene sequences from Athabasca oil 
sands formation waters. For each of the bacterial OTUs the most closely related sequences from GenBank and other related sequences of 
interest (e.g. from other petroleum reservoir habitats) are indicated. Bootstrap values next to branching nodes are based on 100 resamplings. 
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Fig. 5. Concentrations of dibenzothiophenes (the sum of Co-2-alkyl 
dibenzothiophenes) in crude oil (white squares) and of Co-5-all<yl 
naphthalenes in crude oil (black circles) as a function of reservoir 
depth in Athabasca oil sands at the Muskeg River mine. The 
profiles are indicative of microbial degradation of crude oil 
compounds catalysed at the oil water transition zone over 
geological timescales (Head etal., 2003). Formation waters were 
obtained from dewatering wells penetrating the basal water leg at 
c. 80 m depth. 

To assess the possibility that organic sulfur compounds 
provide electron donors for oil-associated Epsilonproteo- 
bacteria (as reported for S. kujiense strain YK-1 ; Kodama 
and Watanabe, 2003) dibenzothiophene concentrations 
were measured at different depths in drill core samples 
obtained from the oil sands layer overlying the formation 
waters. Figure 5 shows the dibenzothiophene concentra- 
tion in oil decreasing from 215 ppm near the surface to 4 
ppm near the oil-water contact at 80 m. This profile 
mirrors the general trends of oil biodegradation with 
depth, as shown by the concentrations of Co to C5 alkyI 
naphthalenes, which decrease from 1011 ppm to 12 ppm 
in the same samples (alkyI aromatic hydrocarbon distri- 
butions and concentrations are indicative of the extent of 
biodegradation in heavy oils; Head etal, 2003; Larter 
etal., 2008). These depth profiles are consistent with 
microbial activity occurring over geological timescales at 
the oil water contact (c. 80 m depth in this system; further 
details provided in Fustic etal., 2011), and suggest that 
microbial desulfurization of crude oil compounds has 
occurred in these oil sands. 



Discussion 

An analysis of 19 clone libraries of bacterial 16S rRNA 
genes from different oil fields revealed that Epsilonpro- 
toeobacteria are among the most frequently detected taxa 
in fossil fuel reservoir production fluids (Fig. 6). On the 
basis of these published reports, only Firmicutes and 
Gammaproteobacteria have been detected in oil fields 
more often than Epsilonproteobacteria. Grouping these 
19 studies according to low-temperature (< 50°C) and 
high-temperature (> 50°C) reservoirs reveals that in low- 
temperature systems like the near-surface oil sands res- 
ervoir investigated here, Epsilonprotoeobacteria are the 
most frequently reported bacterial group (Fig. 6A). Fur- 
thermore, in low-temperature oil fields in western Canada, 
Epsilonprotoeobacteria are represented at high abun- 
dance in 16S rRNA gene clone libraries (Fig. 6B), as 
demonstrated here for the Athabasca oil sands (Table 2). 
By applying multiple bacterial PCR primers to DNA 
extracted using different techniques, our study confirms 
that Epsilonproteobacteria are indigenous to low- 
temperature petroleum reservoirs where they are abun- 
dant members of the microbial community; their massive 
dominance inferred from several 16S rRNA gene-based 
analyses is unlikely to be the result of primer selectivity as 
has been previously suggested (Grabowski etal., 2005a; 
Pham et ai, 2009; Head et al, 201 0). The oil sands at the 
Muskeg River mine are part of the same Athabasca heavy 
oil system as the Pelican Lake reservoir (Fig. 1). The 
recovery of 16S rRNA gene sequences from Epsilonpro- 
teobacteria (99-100% of clones) in pristine formation 
waters from both oil fields is in good agreement, and 
suggests that the dominance of Arcobacter spp. in the 
Pelican Lake formation waters (Grabowski etal, 2005a) 
reflects an important role or roles for Epsilonproteobacte- 
ria in heavy oil reservoirs. 

The most abundant phylotype in all three of our bac- 
terial libraries belonged to the genus Sulfuricurvum. 
Unlike the other epsilonproteobacterial genera detected 
{Arcobacter and Sulfurospirillum spp.; Table 2; Fig. 4), 
Sulfuricurvum spp. are not known to be capable of 
chemoorganotrophic metabolism on the basis of sub- 
strate tests with the only cultivated member of this 
genus, S. kujiense strain YK-1 (Kodama and Watanabe, 
2003; 2004). Despite being autotrophic, strain YK-1 
grows particularly well by oxidizing reduced, most prob- 
ably organic, sulfur compounds in crude oil. This trait 
may explain Sulfuricurvum dominating an oil sands envi- 
ronment with a relatively high sulfur content (Strausz and 
Lown, 2003). The oil-water contact in the Muskeg River 
oil sands shows distinct gradients of dibenzothiophenes 
indicative of biodegradation in the aquifer associated 
with the oil sands (Fig. 5; Head etal., 2003). Biological 
sulfur oxidation in these oil sands is apparent from the 
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Fig. 6. Results from 19 published bacterial clone libraries (references are given in B) from subsurface fossil fuel reservoirs were evaluated for 
the occurrence of dominant taxa. Seven of the habitats were low-temperature environments (reported as < 50°C in situ) as Indicated by the 
dashed line in (A), and 12 were high-temperature oil reservoirs (> 50°C). (A) shows that Epsiionproteobacteria are the third most frequently 
occurring group overall and were detected In five out of seven low-temperature systems. The abundance of Epsiionproteobacteria In these oil 
field clone llbrahes is summarized in (B), which indicates that Epsiionproteobacteria are particularly abundant In low-temperature oil fields In 
western Canada. The hatched bar indicates the discrepancy in results reported from the Kuji oil storage cavity (Watanabe etal., 2000; 2002) 
as described in the introduction. 



■ >50°C 
□ <50°C 



isotoplc composition of sulfate (mean -i-23.3%o) being 
enriclied in ^^S relative to elemental sulfur (mean 
-1-26. 7%o) in the formation waters, even in the absence of 
nitrate (Table 1) - the electron acceptor that supports 
growth of S. kujiense on organic sulfur compounds in oil 
(Kodama and Watanabe, 2003). Oxidation of reduced 
sulfur compounds in situ may have been driven by micro- 
bial reduction of metal species (Aller and Rude, 1988; 
Beal etal., 2009), and the presence of Fe^+ and Mn^+ in 
the formation waters (Table 1) is indicative of metal 
cycling in this system. Many Arcobacter and Sulfurospir- 
illum spp. have been physiologically characterized and 
some are capable of iron or manganese reduction 
(Thamdrup etal., 2000; Fedorovich etal., 2009; Fry 
etal., 2009); however, S. kujiense strain YK-1 has not 



been tested for this capability. Like some other Epsiion- 
proteobacteria, strain YK-1 can grow under microaero- 
philic conditions (Kodama and Watanabe, 2003; 2004), 
and although these oil sands formation waters were 
noticeably sulfidic during sampling, rapid metabolism of 
low levels of oxygen entering this system via groundwa- 
ter flow cannot be excluded. 

Sulfuricurvum spp. have been detected in many 
groundwater environments (Campbell etal., 2006; Porter 
and Engel, 2008; Fig. 4). The Kuji underground crude oil 
storage cavity (the source of strain YK-1; Kodama and 
Watanabe, 2004) is not a pristine oil reservoir like the oil 
sands studied here, but rather an engineered subsurface 
habitat where temporarily stored crude oil comes into 
contact with the local groundwater (Watanabe etal., 
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2000). Indeed, Suifuricurvum spp. have not been reported 
in oil field microbial diversity studies (e.g. those indicated 
in Fig. 6) until now. The subterranean biogeography of 
several Suifuricurvum spp. in groundwaters (Porter and 
Engel, 2008) suggests that a physical connection may 
exist between this Athabasca oil sands reservoir and the 
surrounding basal aquifer, which would influence the oil 
sands reservoir microbiota and biogeochemistry at oil 



water contact zones where biodegradation occurs (Head 
etal., 2003). The ways in which indigenous petroleum 
reservoir microbial communities establish are largely 
unknown, and groundwater infiltration - thought to influ- 
ence microbial community composition in some coal 
deposits (Tseng etal., 1998; Schlegel etal., 2011) - may 
influence the resident microbiota in shallow western 
Canadian biodegraded oil sands. 
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Hydrogenotrophic methanogenesis is thought to be 
responsible for the formation of heavy oil (Jones etal., 
2008), and depends on CO2 and H2 production via fer- 
mentation of crude oil components. It remains unclear 
which taxa catalyse syntrophic fermentations in biode- 
graded subsurface petroleum reservoir environments 
(Pham etal., 2009), and the predominance of Epsilonpro- 
teobacteria is intriguing in this context. Certain studies 
suggest Epsilonproteobacteria might be involved in 
syntrophic anaerobic communities that degrade hydrocar- 
bons. For example, an epsilonproteobacterium was 
recently shown to assimilate ^^C during syntrophic degra- 
dation of ^^C-labelled benzene in sulfate-reducing enrich- 
ment cultures (Herrmann etal., 2010) and appears to 
contribute to benzene fermentation, e.g. into H2, acetate 
and CO2. Indeed in the study of Herrmann and colleagues 
(2010), the highest "C enrichment was observed for the 
Epsilonproteobacteria. Sulfurospirillum spp. can ferment 
a wide variety of organic compounds (Luijten et al., 2003) 
and were among the fermentative heterotrophs isolated 
from a subsurface coal deposit (Fry etal., 2009) and the 
Pelican Lake oil reservoir (Grabowski etal., 2005a; 
Fig. 4). Acetate is metabolized by Arcobacter spp. in 
various anoxic settings (Thamdrup etal., 2000; Fedorov- 
ich etal., 2009; Webster etal., 2010). Syntrophic oxida- 
tion of acetate to H2 and CO2 is important in some 
methanogenic crude oil degrading systems (Jones etal., 
2008), and is catalysed by a few cultured isolates (Hattori, 
2008; Westerholm etal., 2010) that can also catalyse the 
reverse reaction (i.e. acetogenesis from H2 and CO2; Lee 
and Zinder, 1988). Several Epsilonproteobacteria can use 
H2 as an electron donor (e.g. Gevertz etal., 2000; 
Kodama etal., 2007), including the acetogenic /4rco- 
bacfer dominating the bacterial community in the Pelican 
Lake reservoir (Grabowski etal., 2005a). This is consis- 
tent with Arcobacter sequences dominating the metage- 
nome of formation water from a coalbed methane 
reservoir (US Patent 2010/0047793 Al, 2010), a habitat 
where syntrophic methanogenic hydrocarbon degradation 
also occurs (Strgpoc etal., 2008). 

Intact polar lipid analyses indicated that bacteria pre- 
dominated over archaea in the formation waters (Olden- 
burg etal., 2009). Taken together, lipid and DNA-based 
results therefore suggest a Sulfuricurvum sp. is the 
most abundant organism in this heavily biodegraded oil 
reservoir, with Arcobacter, Sulfurospirillum and hydro- 
genotrophic Methanomlcroblaceae spp. making up 
smaller fractions of the microbial community (Table 2). 
The latter groups may be involved in methanogenic crude 
oil degradation that has generated these and other heavy 
oils, and currently these populations may have given way 
to a dibenzothiophene-degrading Sulfuricurvum popula- 
tion well suited to this sulfur-enriched oil sands reservoir 
(around 5% sulfur by weight; IVIarcano, 2011). Our results 



provide clear evidence that 16S rRNA gene-based detec- 
tion of Epsilonproteobacteria in oil field formation waters 
is not artefactual. On the contrary, Epsilonproteobacteria 
are predominant members of low-temperature biode- 
graded petroleum reservoirs, and their ecophysiology in 
such environments can be explained by examining differ- 
ent aspects of crude oil sulfur and hydrocarbon bio- 
geochemistry. 

Experimental procedures 

Sample collection 

Athabasca oil sands formation waters were sampled at the 
Muskeg River oil sands mine north of Fort McMurray, Alberta, 
Canada (Fig. 1). Mining operations at this site require reser- 
voir formation water discharge in advance of open pit oil 
sands removal and water from three out of the six dewatering 
wells surrounding the active excavation area were sampled 
(Fig. 1). Wellhead valves were opened and flushed for 
several minutes prior to collecting 1 2 or 24 I using a sterilized 
metal funnel to completely fill sterile 4 I brown glass bottles 
that had previously been flushed with nitrogen gas. Water 
samples were transported to the laboratory in a cool box 
where they were subsequently stored at 4°C. Oil sands 
samples were obtained from areas of the reservoir that had 
been undergoing commercial excavation. Large volumes 
(several m^) of sediment were excavated from the intact 
formation by mine equipment (large backhoe). This enabled 
aseptic subsampling from the middle part of the excavated oil 
sands. Subsamples were stored anaerobically in N2-flushed 
containers kept at room temperature. Oil sands were also 
obtained from vertical drill cores spanning the oil leg and 
reaching depths close to the oil water transition zone. 



Geochemical analysis 

Formation water samples from dewatering wells numbered 1 , 
4 and 6 were subject to a suite of geochemical analyses, 
performed by the Analytical Geochemistry Group at the Uni- 
versity of Calgary. To determine isotopic compositions of 
sulfate and elemental sulfur in the formation waters, BaS04 
was precipitated following addition of 0.5 M BaCl2 and 
elemental sulfur was extracted following filtration and Bligh 
and Dyer extraction (Bligh and Dyer, 1959). Sulfur isotope 
ratios were determined by continuous flow-elemental 
analysis-isotope ratio mass spectrometry (CF-EA-IRMS) with 
a Finnigan Mat Delta+XL spectrometer interfaced with a 
Carlo Erba NA 1 500 elemental analyser and are expressed 
relative to the international standard Canyon Diablo Troilite. 
For IPL analysis, biomass was concentrated by vacuum fil- 
tration of the 12 or 24 I samples through 0.2 |xm pore size 
filters (Millipore) that were subsequently freeze dried to allow 
biomass removal by scraping it from the filters with a steril- 
ized metal spatula. The biomass was then subjected to a 
modified Bligh and Dyer extraction (Bligh and Dyer, 1959) to 
recover polar fractions which were analysed by liquid 
chromatography-mass spectrometry to determine IPL com- 
position, as described previously (Oldenburg etal., 2009). 
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Bitumen samples from the oil sands drill core were sub- 
jected to petroleum geochemistry analyses. 'Aromatic hydro- 
carbon' fractions including dibenzothiophenes were isolated 
from de-asphaltened bitumen extracts following the proce- 
dure of Bennett and Larter (2000). AlkyI naphthalenes were 
quantified relative to a D8-naphthalene standard (response 
factors of unity were employed therefore data can be consid- 
ered semi-quantitative). Mass spectral characterization of 
compounds in the aromatic hydrocarbon fractions was 
carried out using splitless injection GC MS with an Agilent 
6890 chromatograph interfaced to a 5973 quadrupole mass- 
selective detector. Saturated and aromatic hydrocarbons 
were analysed with a DB-5-coated fused silica column (30 m 
length, 0.32 mm id, 0.25 |im film thickness; J&W Scientific) 
according to the following temperature programme; 40°C 
(held for 2 min) to 300°C (held for 20 min) increasing at 
4°C min ^ 



DNA extraction and PCR amplification of 16S 
rRNA genes 

Most of the blomass Isolated by filtration was dedicated to I PL 
analysis as described above. DNA was extracted from the 
remaining blomass residue on the same filters, which allowed 
for 16S rRNA gene analyses. Initial screening was performed 
by DGGE to confirm that the microbial communities in the 
three water samples were similar. Subsequent community 
analyses were performed on water from wellhead #1 accord- 
ing to two complementary approaches: 

(i) At the University of Calgary, filters were submerged in a 
resuspension buffer (0.15 M NaCI; 0.1 M EDTA) and 
Incubated on a rotary shaker overnight. DNA was 
extracted from resuspended blomass using the method 
of Marmur (Marmur, 1961) with modifications (Voordouw 
etal., 1990) including successive rounds of freezing and 
thawing as well as lysozyme, RNase and proteinase K 
treatments. DNA concentrations were measured by 
spotting dilutions on a square of 1% agarose gel con- 
taining ethidium bromide (c. 10|a,g mT'). The concentra- 
tion of DNA was estimated under UV illumination by 
comparison with fluorescence of known amounts of 
phage lambda DNA. The presence of bacterial DNA was 
confirmed in all formation water samples by positive 
PGR amplification from all DNA extracts using primers 
8f and 1406r (see Table SI). Extracted DNA from well- 
head #1 formation water was sent to the JCVI in Rock- 
ville for amplification of 16S rRNA genes from Archaea 
using primers arch8f and arch1492r, and from Bacteria 
using primers 9f and 1545r (Table SI). Amplifications 
were performed using a DNA Engine Tetrad PTC-225 
thermal cycler (MJ Research, Waltham, MA) with an 
Initial denaturation of 2 min at 94°C, followed by 29 
cycles of 30 s at 94°C, 30 s at 55°C and 2 min at 72°C, 
with a final extension of 5 min at 72°C. A negative 
control PCR reaction in which the genomic DNA tem- 
plate was replaced by an equivalent volume of sterile 
distilled water was also included. 

(II) Parallel analysis was conducted at Newcastle University 
where DNA was extracted directly from the filters for each 
wellhead formation water sample using the FastDNA Spin 



Kit for soil (Q-BIOgene, UK). Filters containing blomass 
residue were cut into small sections and added to 2 ml 
tubes and DNA was extracted following the manufactur- 
er's guidelines, which include an initial bead-beating step 
using a RiboLyzer (Hybaid). 16S rRNA genes from 
Archaea were amplified by PCR using primers arch46f 
and arch1017r. 16S rRNA genes from Bacteria were 
amplified using two different primer pairs, either 8f and 
1542r or inosine-341f and 1492r (Table 2). All primer 
sequences and a summary of downstream procedures 
used to analyse the amplified 16S rRNA genes are listed 
in Table SI. 



Denaturing gradient gel electrophoresis (DGGE) 

Shorter 16S rRNA gene fragments were amplified for DGGE. 
For bacterial 16S rRNA genes a PCR product generated 
using primers 341f-GC and 534r was used directly for DGGE. 
Bacterial DGGE was also performed following nested PCR 
amplifications where PCR products generated with inosine- 
341 f and 1492r or 8f and 1542r were used as templates for a 
second PCR reaction with 341f-GC and 534r. A nested 
approach was also used to analyse archaeal 1 6S rRNA gene 
fragments by DGGE. PCR products generated with primers 
arch46f and arch1017r were used as template for a second 
round of PCR using primers arch344f-GC and Uni522r. All 
primer sequences are listed in Table SI. DGGE was con- 
ducted in 10% acrylamide gels with 30-60% denaturant, as 
described previously (Rowan etal., 2003). Acrylamide gels 
were run for 4 h at 200 volts using the D gene system (Bio- 
Rad, Hemel Hempstead, UK) and were subsequently stained 
for 30 min in SYBR green I (Sigma, Poole, UK; diluted 
1/10000 in 1x TAE). Stained gels were viewed under ultra- 
violet light and gel images were recorded using a Bio-Rad 
Fluor-S® Multilmager (Bio-Rad, UK). Bands of interest were 
excised for sequencing. 

Cloning, sequencing and phylogenetic analyses 

Bacterial and archaeal 16S rRNA gene fragments were 
cloned using TOPO TA cloning kits (Invitrogen) according to 
the manufacturer's instructions. Approximately 384 clones 
were analysed from the JCVI libraries, and approximately 96 
clones were analysed from each of the Newcastle libraries. 
Nearly complete full-length 16S rRNA gene sequences 
{Escherichia coli positions 9-1545) were obtained from the 
JCVI libraries by sequencing using Ml 3 cloning vector 
primers. Partial 16S rRNA gene sequences were obtained 
from Newcastle clones using internal 16S rRNA gene 
primers (Edwards etal., 1989). Closely related sequences 
from GenBank were identified by blast searching (Altschul 
etal., 1990) of the GenBank database and using the Ribo- 
somal Database Project SeqMatch tool (Cole etal., 2007). 
Sequences were aligned using the SILVA web aligner 
(Pruesse etal., 2007) and alignments were manually cor- 
rected using BioEdit (Hall, 1999). Neighbour joining phylo- 
genetic trees were constructed using MEGA4 (Tamura etal., 
2007). Sequence alignments for phylogenetic tree recon- 
struction included the top three hits obtained in blast and 
SeqMatch searches (multiple environmental sequences 
from the same study or habitat were considered as one hit) 



©2011 Society for Applied Microbiology and Blackwell Publishing Ltd, Environmental Microbiology, 14, 387-404 



400 C. R. J. Hubert et a\. 



to each of the oil sands phylotypes, as well as other 
sequences of interest. Representative sequences from the 
formation water clone libraries have been deposited in the 
GenBank database under Accession No. JF789587 to 
JF789598. 

Meta-analysis of published oil field 16S rRNA bacterial 
clone library results 

A literature survey revealed 19 published clone libraries of 
bacterial 16S rRNA genes from oil reservoir production fluids. 
Clone library results were collated to reveal the proportion of 
clones affiliated with major bacterial divisions. This analysis 
considered over 3600 sequences of cloned 1 6S rRNA genes, 
which we classified at the level of major phylogenetic groups 
(phylum/subphylum level). From these data, distributions 
based on the frequency of occurrence and average percent- 
age representation of major taxa in clone libraries could be 
plotted (see Fig. 6). Only clone libraries that used broad 
specificity ('universal') primers for PGR amplification of 16S 
rRNA genes were included in the survey. Published DGGE 
analyses on oil field production fluid samples were not 
included in this analysis. Results were considered in the 
context of low- and high-temperature petroleum reservoirs 
using a cut-off of 50°G {in situ reservoir temperature as indi- 
cated in the publications). We recently used a similar 
approach to conduct a broad assessment of microbial com- 
munity structure in petroleum-impacted environments includ- 
ing soils and surface sediments (Gray etal., 2010). 
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